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Abstract 



An explicit formula for the ADM mass of an asymptotically AdS black hole in a 
generic Lovelock gravity theory is presented, identical in form to that in Einstein 
gravity, but multiplied by a function of the Lovelock coupling constants and the 
AdS curvature radius. A Gauss' law type formula relates the mass, which is an 
integral at infinity, to an expression depending instead on the horizon radius. 
This and other thermodynamic quantities, such as the free energy, are then 
analyzed in the limits of small and large horizon radius, yielding results that are 
independent of the detailed choice of Lovelock couplings. In even dimensions, 
the temperature diverges in both limits, implying the existence of a minimum 
temperature for black holes. The negative free energy of sufficiently large black 
holes implies the existence of a Hawking-Page transition. In odd dimensions 
the temperature still diverges for large black holes, which again have negative 
free energy. However, the temperature vanishes as the horizon radius tends to 
zero and sufficiently small black holes have positive specific heat. 



1 Introduction 



Many interesting, stationary black holes are not known in analytic form. The list includes 
charged rotating black holes in D > 4, black rings in D > 5, localized Kaluza-Klein black 
holes and rotating black holes in Lovelock gravity theories. Results have been obtained 
using a variety of approximate techniques, including e.g. perturbative expansions in the 
slowly rotating limit [U El El HI E] or effective field theory methods [SI El E]- However, one 
may also ask whether in the absence of analytic solutions, it might still be possible to at 
least obtain certain properties of such spacetimes, such as thermodynamic ones, exactly. 

In this paper we will address this question for a class of spacetimes that might best be called 
"semi-unknown" , namely static Lovelock black holes. It is known [9] that these spacetimes 
are Schwarzschild-like, in the sense that the metric written in the general static, spherically 
symmetric form 

dv^ 1 

ds 2 = -4>(r)dt 2 + J^ + r 2 dQ 2 D _ 2 , (1) 

additionally has <j>(r) = f(r). For a generic Lovelock theory, the function f(r) must satisfy 
an algebraic equation of order [=^1 , where D is the spacetime dimension and the closed 
brackets denote taking the integer part. The coefficients of the polynomial are the coupling 
constants of the higher curvature terms in the Lovelock Lagrangian. Since the roots of a 
generic polynomial equation can be found in terms of radicals only up to order 4, it follows 
that in spacetimes dimensions D > 10 the general solution for a static Lovelock black hole 
cannot be written down in a simple closed form and hence remains unknown^. 

We make use of results from references [T8 ] [T9 | [20]. derived using the Hamiltonian formula- 
tion [21] . to obtain exact properties of Lovelock black holes, without requiring the explicit 
(and in general unknown) form of the metric function /(r). In fact with future applica- 
tions to other, even more unknown, black holes in mind, we will "forget" that we know 
that 0(r) = f(r) and take the general static, spherically symmetric form ([1]) as our starting 
point. Our focus throughout will be on black holes with asymptotically AdS boundary 
conditions. 



Further motivation comes from recent work on the CFT duals of asymptotically AdS so- 
lutions to higher curvature gravities. Significant work has been done, for example, on the 
relation between CFT plasmas and their gravitational duals, [221 1231 12H [25] and on im- 
plications for the CFT of causality and stability in the bulk using holography [261 ETl I2H1 
[291 EQl EH E21 [331 EH [351 ESI EZl EH [39l SO]. As the correspondence between the higher 

1 Given that the solutions to generic cubic and quartic equations are quite cumbersome, in practice the 
general solution is only known in useful form for D < 6, although certain non- generic solutions such as those 
for Gauss-Bonnet gravity [TU1 [TTJ [T2J [T3] , or Chern-Simons gravity [TU [T3] are known in all dimensions. 
Note also that an analytic solution for a rotating black hole in D = 5 Gauss-Bonnet gravity has also been 
found in the case of Chern-Simons couplings [16] (see also |17|). 
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curvature bulk theory and the boundary CFT becomes better understood, it is interesting 
to ask whether the constructs we will make use of below, the Killing-Lovelock potentials 
and the associated Gauss' law relations that connect behavior at the horizon to behavior 
at infinity, have an analogue in the dual CFT. 

The paper will proceed as follows. In section ([2]) we recall the basic elements of Lovelock 
gravity theories [U] that will subsequently be used. In section (j3J) we will derive the 
expression for the ADM mass of an asymptotically AdS Lovelock black hole in terms of 
its far field behavior. The ADM boundary integral receives contributions from each higher 
curvature term as well as from the Einstein term. The formula for the ADM mass is implicit 
in the work of |42j which computed the entropy in Lovelock theories and established the 
first law. However significant additional steps are required to make the result of jl2] fully 
explicit. The expression we find ultimately takes the simple form of the standard ADM 
mass integral multiplied by a function of the Lovelock couplings and the asymptotic AdS 
curvature radius. This formula for the mass constitutes a new result. 

In section we review certain elements from our previous work rj8[[T9j[20] that we will be 
making use of. This includes the Killing-Lovelock potentials which allow for the derivation 
of a Gauss' law type expression for the Hamiltonian constraint, which will be our main 
tool. This expression relates certain boundary integrals at infinity to surface integrals on 
the horizon, without requiring the explicit form of the metric in between. For reference, 
we also present the Smarr formula [20] which gives the mass in terms of the entropy and 
additional thermodynamic quantities that arise in the Lovelock theory. 

The bulk of our new results appear in sections (J5j) and (Q. In section (jSJ) we use the Gauss' 
law formula to obtain an expression for the mass, originally given by the far field behavior 
of the metric, in terms of the horizon radius of the black hole. With the goal in mind of 
providing a simple expression for the free energy, we also present formulas for the entropy, 
surface gravity and other thermodynamic quantities in terms of the horizon radius. The 
detailed behavior of these expressions over the full range of horizon radii depends on the 
choice of Lovelock coupling constants. However, we find that they have generic behavior in 
the limits of small and large horizon radii. There are qualitative differences between even 
and odd dimensions that arise in the small black hole regime. In particular, one finds that 
in odd dimensions there exists a "mass gap", in the sense that the mass tends to a finite 
value as the horizon radius tends to zero. The Hawking temperature also vanishes in this 
limit and sufficiently small black holes have positive specific heat. For even dimensions, 
however, the mass tends to zero for vanishing horizon radius and the specific heat for small 
black holes is negative. 

In section ([6]) we assemble ingredients to construct a relatively simple expression for the 
free energy of black holes in generic Lovelock theories. The free energy had previously been 
computed only up to inclusion of the quadratic Gauss-Bonnet term (T2J [13] in the Lovelock 
Lagrangian. The free energy for a stationary black hole solution is generally obtained 
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by computing its Euclidean action. We emphasize that our method does not require the 
explicit functional form of the metric. We analyze the free energy and also the specific 
heat in the small and large black hole limits and comment upon the Hawking-Page phase 
transition for generic AdS-Lovelock black holes. In section ([7]) we briefly summarize our 
results and offer some directions for further related work. 



2 Lovelock gravity 



The Lagrangian of a Lovelock gravity theory in D spacetime dimensions is given by 

k=0 

where k max = [(D — l)/2] and the are real- valued coupling constants. The symbol 
stands for the contraction of k powers of the Riemann tensor given by 

r(k) _ ]_caib 1 ...a k b k p cidx p c k d k /q\ 
— 2 fc Cldl "-- Cfcdfe ■■■ jrL a k b h ■ {.<>) 

where the 5-symbol is the totally anti-symmetrized product normalized so that it takes 
nonzero values ±1. The term gives the cosmo logical constant term in the Lagrangian, 
while £W gives the Einstein-Hilbert term and the quadratic Gauss-Bonnet term. The 
upper bound in the sum (T5]) comes about because vanishes identically for D < 2k and 
turns out to make no contribution to the equations of motion in D = 2k. 

The equations of motion for Lovelock gravity can be written as Q a i > = where 

q\ = bkg{k)ab - ( 4 ) 

k=0 

We will not need the explicit expressions for these quantities. However, it is crucial in what 
follows that each of the tensors in this sum satisfies a conservation law V a (? < - fc - )a b = 0. 

Depending on the values of the coupling constants bk the theory may have anywhere from 
zero up to k ma x distinct constant curvature vacuum solutions. Because we will be focusing 
on asymptotically AdS black holes below, our considerations will implicitly be limited to 
the subset of theories admitting at least one constant negative curvature vacuum. We will 
denote the curvature radius of this AdS vacuum by /. 

We will also need the Hamiltonian formulation of Lovelock gravity which was developed in 
|21j . As usual in the Hamiltonian picture the spacetime metric is split according to 

9ab = -n a n b + s ab (5) 
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where n a is the unit timelike normal to a spatial slice £ with induced metric s a b and 
these satisfy the orthogonality relation s a t,n b = 0. As in Einstein gravity the time-time 
and time-space components of the field equations act as constraints on initial data. In 
Lovelock gravity the Hamiltonian and momentum constraint operators H = — 2n a n b Q a b 
and H a = —2s a b n c Q^ c are given by the sums 

kmax kmax 

H=Y,hHW, H a =J2hH^ (6) 

k=0 k=0 

with H ^ = —2n a n b Q^j and Ha = —2s a b n c ■ We will particularly need the expression 
for 

Tj(k) _ _ ^ caibi...a k b k p Cl d± p c k d k 

~ ofc Cldl — C * <i * 161 ■■■ lx akbk \' J 

where the tilde on the 5-symbol indicates that its indices are projected with the spatial 
metric sV The actual Hamiltonian function for evolution of initial data with respect to 
a vector field £ a is then given by = F H + F a H a , where the lapse and shift (F,F a ) 
are the components of the vector field £ a normal to and along the spatial slice, so that 
£ a = Fn a + F a . As in Einstein gravity the lapse and shift are Lagrange multipliers. 

Finally, for use below we define the sums 

_1 }k+n 

(k-n)\ l 2k 



v» = E(-D t+ "7F^TW < 8 > 



k=0 



where the rescaled coefficients i>k are given by i>k = (D — 1)! bk/(D — 2k — 1)!. Note that 
the combinations ^/Z 2 ^ -1 -* are dimensionless. The condition S(o) = 0, in fact, determines 
the allowed constant curvature vacua of the theory, while the sums S(i) and sp) will turn 
up in our results below. 



3 ADM Mass 



An expression for the ADM mass of an asymptotically AdS Lovelock black hole can be 
obtained using the methods of Regge and Teitelboim [H]. We find a simple, explicit 
formula which is similar to that for the ADM mass. As in Einstein gravity, a boundary term 
must be included in the Lovelock Hamiltonian to ensure that the Hamiltonian variational 
principle correctly yields the equations of motion. The variation of this boundary term 
cancels another boundary term arising via integration by parts from the variation of the 
bulk Hamiltonian. The ADM mass is defined to be the value of the Hamiltonian divided 
by a factor of lQnG. Since the bulk Hamiltonian vanishes on solutions, the ADM mass is 
simply proportional to the value of the Hamiltonian boundary term. Implicit here is that 
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the vector £ a determining the direction of Hamiltonian evolution should be asymptotic to 
the time translation Killing vector of AdS. 



This procedure was formally carried out in reference [12] as part of establishing the first 
law for Lovelock black holes. Each higher curvature term in the Lovelock theory makes a 
contribution to the Hamiltonian boundary term. For asymptotically flat solutions, possible 
if b = in the Lagrangian ([2"]), because of the fall-off of the curvature tensor, only the 
boundary term corresponding to the Einstein term is nonzero at infinity. The formula for 
the mass is then the same as in Einstein gravity. However, the non-zero asymptotic value 
of the curvature tensor for AdS Lovelock black holes leads to contributions from all the 
higher curvature boundary terms. We will see, in fact, that the result is given by the usual 
ADM integral multiplied by a function of the Lovelock couplings and the background AdS 
curvature. To our knowledge this expression has not been derived in the literature^). 

The Regge-Teitelboim type construction of the ADM mass in [12] yields an expression of 
the form 



We will determine the boundary integrands B^ c explicitly for asymptotically AdS bound- 
ary conditions. Because our interest is focused on static black holes, we will assume that 
the momentum and the shift vector vanish sufficiently rapidly in the asymptotic region that 
they do not play a role in this construction. It is then sufficient to retain only the term 
"H^ = FH in the Hamiltonian and to take the curvature tensor in (J7J) to be that of the 
spatial metric s a b. Now, assume that the spatial metric s a b solves the Lovelock constraint 
equations with asymptotically AdS boundary conditions and add to it an arbitrary pertur- 
bation h a t, that also respects these boundary conditions. To first order, the perturbation 
to the Riemann tensor is then SR ab cd = R ab hf* — 2D[ a D^°h b , d \ where D a is the covariant 
derivative operator for the spatial metric s ab . Plugging this in to the variation of (J7]) yields 
the perturbation of the feth Hamiltonian 



= -7T k S.-.-.tl Ra lbl ^ ■ ■ ■ Ra k _ lbk -i Ck - ldk - 1 {Ra^he^ ~ 2D a . k D^h b ^) (10) 



2 A number of related results have been established in the literature. The Hamiltonian approach of 
[44] has been used to obtain an expression for the mass in pure Lovelock gravity theories [15] with only 
a single term in Lagrangian. Significant work has also been done in the case of Gauss-Bonnet gravity. 
Deser and Tekin have extended the formalism of Abbott and Deser [35] to general quadratic theories of 
gravity [46, 47]. The mass in Gauss-Bonnet gravity has also been found by means of a Noether's current 
construction by Deruelle et. al. [48] and using the Palatini formalism by Katz and Livshits [49] The mass 
of asymptotically AdS black holes in general higher curvature theories has also been addressed in [3JJ] • 

Padilla [ST] uses the Cartan formalism to derive a general expression for the mass in Gauss-Bonnet 
gravity by transforming the trace K boundary term in the action into Hamiltonian variables. By assuming 
that a generalization of this result holds in general Lovelock theories, the mass of static Lovelock black 
hole is identified by Cai [52] with a constant of integration arising from integration of the Hamiltonian 
constraint. Our general derivation of the mass in Lovelock theories demonstrates the correctness of this 
assumption. 
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Combining these, multiplying by the lapse function F, integrating by parts, and making 
use of the Bianchi identity for the Riemann tensor then gives an overall expression of the 
form 

F6H = 6H* ■ F + ^b k D c B {k)c . (11) 

k 

Here the first term represents an adjoint differential operator acting on the lapse function 
F and the vectors B^ c are given by 

B [k)c = |- Rm ^ eih ■ ■ ■ R ^^ ek - lfk - 1 ( FDah « ~ h * bDaF ) ( 12 ) 

These are the vectors appearing in the expression for the ADM mass given above. The 
quantity h a b is then the deviation of the static black hole metric from the asymptotic AdS 
background, while the minus sign in (jUJ) arises from the cancellation of boundary terms 
required in the Regge-Teitelboim construction. The expression can be made more fully 
explicit by noting that near infinity only the leading order background AdS curvature 
R a b cd = (V^ 2 )^afe contributes. This leads to the result 

M = (o-ik'o-2) ' (isb L da < - m ' b) - hD ' F + h ' bD ^) ' (13) 

The part of the formula in parenthesis is the usual expression for the ADM mass. The sum 
S(i), defined in (J5J), depends on the Lovelock coupling constants and the curvature radius 
I of the asymptotic AdS vacuum. For Einstein gravity with a cosmological constant the 
prefactor outside the parenthesis reduces to one. The asymptotic form of the lapse function 
is F = yl + r 2 /l 2 , and therefore terms in (1131) involving the derivative of the lapse can 
make non-trivial contributions to the mass. In the asymptotically flat case, these terms 
fall off too quickly to contribute and the integrand reduces to that of the ADM mass in the 
asymptotically flat case. 



4 Killing-Lovelock Potentials 



Killing-Lovelock potentials [T8| [20] allow the Lovelock Hamiltonian constraint equations 
for a spacetime with a Killing vector to be written in a Gauss' law form. In reference [20J 
this property was used to derive an extended first law including variation of the Lovelock 
coupling constants and also a related Smarr formula for Lovelock black holes. This will 
also be our primary tool below. 

For a spacetime, such as a stationary black hole, with a Killing vector £ a the Killing- 
Lovelock potentials were defined in [18] to be antisymmetric tensors (3^ ab satisfying the 
relations 

Va/3 (k)ab = _ 2 g(Q* b £b (14) 
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Their existence is guaranteed by the vanishing divergence of the right hand side. However, 
they are not uniquely determined, since a divergenceless tensor such as the divergence of 
an arbitrary 3-index antisymmetric tensor may always be added. This ambiguity does not 
affect results such as the extended first law or Smarr formula. 

Now consider Hamiltonian evolution with respect to the Killing vector, so that the lapse and 
shift are given by the decomposition £ a = Fn a + F a . The definition of the Killing-Lovelock 
potentials allows us to write the kth Hamiltonian function as a total divergence, 

U (k) = FH (k) + F a R {k) = _ 2 g{k)d^c nd = Dc ^{k)cd nd y ^ 

The Hamiltonian constraint equation T-t — can then be written in the form of Gauss's 
law. It ultimately follows from this (see [20] for how this works) that the Killing-Lovelock 
potentials make contributions to the thermodynamics of AdS-Lovelock black holes that are 
similar in form to that of the electrostatic potential in the case of a charged black hole. 
For example, one contribution to the change in energy in the extended first law of [20] 
is proportional to the difference in the integrals of Killing potential over the boundaries 
of a spatial slice at infinity and at the black hole horizon. The AdS boundary conditions 
require a subtraction at infinity of the Killing-Lovelock potential fills' fo r the asymptotic 
AdS spacetime resulting in thermodynamic potentials 

e (k) = - (| dar c (P^ cd - $lf)n d - J dar c ^ cd n d ^ . (16) 

that multiply variations Sb^ of the Lovelock couplings in the extended first law. As a 
consequence of the constraint equations these thermodynamic potentials satisfy the sum 
rule 

£ b k e {k) = o. (17) 

fc=0 

The Smarr formula relates the mass of a black hole which is defined in terms of the behavior 
of the metric near infinity to the entropy which comes from the behavior at the horizon. 
For static Lovelock black holes the Smarr formula [20J, which we will also make use of 
below, is given by 

((D-3) + 2^]m = — [(D- 2)S-S'] -0 (18) 
V 3(1)7 27r 

The entropy S [42J has contributions from the higher curvature Lovelock terms and is given 
by S = A/AG with A = J2 k hA k and 

A k = k [ d^x^C^i^) (19) 
Jh 
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where 7 tt & is the induced metric on the boundary of the spatial slice at the horizon. There 
is also an additional entropy-like contribution at the horizon as well as the net contribution 
Lovelock thermodynamic potentials which are given respectively by 

fc=0 fc=0 

The Smarr formula is derived in [20] via an overall scaling argument from an extended 
form of the first law in which variations of the dimensionful Lovelock coupling constants 
are taken into account. The contribution of S' to the Smarr formula arises from the explicit 
dependence of the entropy on the Lovelock couplings, and the second term in the parenthesis 
on the left in (lisp similarly arises from the dependence of the mass on the couplings b^. 

For Einstein gravity with vanishing cosmological constant one finds that s^) — S' = = 
and one recovers Smarr formula for static, asymptotically flat black holes. The case of 
Einstein gravity with non-zero cosmological constant was discussed in [19]. In this case, 
one still has S(2) = S' = 0. However, there is a non-trivial contribution to from the 
product! A0W. In [T3] it was argued that 0^°^ can be regarded as minus an effective 
volume behind the black hole horizon. Since the cosmological constant is proportional to 
minus the background pressure, one sees that the k = term in the Smarr formula has 
the form pV which is familiar from classical thermodynamics. The mass M should then be 
regarded as the spacetime enthalpy. This interpretation was explored further in [53| [Mj. 



5 Lovelock black holes 



In this section we will derive formulas for the basic thermodynamic properties of static, 
spherically symmetric AdS Lovelock black holes in terms of the horizon radius and 
the Lovelock parameters bf.. We consider static, spherical symmetric, asymptotically AdS 
spacetimes with AdS curvature radius I. The metric can then be taken to have the form 
(PQ), where at large radius the metric functions have the asymptotic forms 

2 2 

^(O-i+p-^s, /(0 ~ i + p - ^ (21) 

for some constants q and c r . We assume that there is a Killing horizon at r = r^, where 
<K r h) = 0- Static spherically symmetric black hole solutions of Lovelock gravity theories 
have been known for some time, starting with the work of PHUT] in the Gauss-Bonnet case 
and [9] in the general Lovelock case (see also [55J). These solutions all have 0(r) = /(r) 
and the field equations reduce to the requirement that /(r) solve a certain polynomial 

3 Onc has A = — 26q in this case. 



8 



equation of order k maX} with coefficients determined by the Lovelock coupling constants^. 
Except in certain special cases, such as Einstein or Gauss-Bonnet gravity or when the 
Lovelock couplings are tuned such that the polynomial has a unique degenerate root [To] , 
the solutions for f(r) are generally not known explicitly 

Here we will use the general relations given in the previous sections to derive expressions 
for the mass M, surface gravity k, entropy S and free energy F of Lovelock black holes 
in terms of without requiring the explicit solution of the field equation. This approach 
provides a geometrical understanding of the formulae. As noted in the introduction, with 
the application to other even less explicitly known black hole spacetimes in mind, we will 
take the general form of the metric ([1]) as our starting point and "forget" that we know 
that solutions to the field equations will have <f)(r) = f(r). 



5.1 Constant curvature vacua and Killing potentials 

Working in the Hamiltonian picture, we assume that for a spherically symmetric static 
black hole the spatial metric has the form 

dv 2 

s ab dx a dx b = -^-y + r 2 dtt 2 D _ 2 , (22) 

that the extrinsic curvature K ab vanishes, and that the Hamiltonian evolution is carried 
out along the static Killing field £ a = Fn a . The function F 2 must then have the same large 
r fall off conditions as the metric function <fr in (|2ip and satisfy F 2 (rh) = 0. Substituting 
the spatial metric into the Lovelock Hamiltonian functions H^ h ' in (I7j) we find that we find 
for these metrics that 

» W = -1^^ (^-^"(l - /)') (23) 
where 7^ = p^-i)! • ^ e constraint equation is 

H=J2 h kH (k) = (24) 

k=0 

First consider the background AdS metric with radius of curvature I and no black hole. 
The AdS functions are F\ dS = f^ds = 1 + r 2 /l 2 . Since the metric satisfies the constraint 
equation with fAds, substituting into ( 124]) gives a relation between the couplings and I 2 , 

k=0 



4 Specifically, if one defines the function F according to / = 1 — r 2 F, then F must satisfy the polynomial 
equation X)fc=(T bkF k — r S-\ for some constant u>. Constant curvature vacua of Lovelock theories solve 
this equation with lj = 0. 
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As noted above, this is the condition that the sum S(o) defined above in (jSJ) vanish in 
a constant curvature vacuum. Whether, or not, there are real positive solutions for I 2 
depends on the values of the Lovelock couplings. We will assume that we are working with 
a Lovelock gravity theory that has at least one real, positive solution I 2 . 

Next we find the projected Killing-Lovelock potentials (3^ rd nd, which are solutions to 
equation ffl5l) with vanishing shift vector F a . Plugging in the form of the kth Hamiltonian 
f[25j) this becomes 

vj£ (^rP (k)rdn «) = - F ^ r (^-^a - m (26) 

Now, recall that the existence of the Killing-Lovelock potentials depends on the the exis- 
tence of a Killing vector £ a , but not on the metric being a solution to the field equations. 
The statement that the thermodynamic potentials O^ sum to zero in ( 1T7|) . however, does 
assume that the constraint equation is satisfied. Since we will want to make use of (ITT)) 
our current task is to solve for the Killing-Lovelock potentials assuming that / solves the 
constraint equation 021]). Note also that the differential equation ( 12B)) includes a factor of 
the lapse function F. In the Hamiltonian picture the lapse is a gauge choice that may be 
freely specified^. We can choose any function F consistent with — F 2 being the norm of the 
static Killing field for a spherically symmetric AdS black hole. Specifically, we can choose 
any F that depends only on r with F(r^) = and satisfying the large r fall off conditions 
in (H]). 

On the other hand, inspection shows that equation ( 126)) is easy to solve if we can choose 
F = y/f. This choice has the correct fall off conditions at large r, since the function f(r) 
does. Therefore, we only need to worry about the condition that F vanishes at the horizon, 
or equivalently whether the function f(r) necessarily vanishes at the horizon? Of course 
it is well known that the Schwarzchild solution in Einstein gravity has F = y/f, with 
f(r) vanishing at the horizon. As noted above the static Lovelock black holes continue 
to have this Schwarzchild form [9]. However, one can also demonstrate that /(r^) = 
under weaker conditions than requiring that the full set of field equations be solved. Let 
us assume that the surface gravity for the black hole metric ([I]) is finite and non-zero, and 
that the constraint equation is satisfied at the horizon radius r^. We show in the Appendix 
that these conditions imply that / ~ F 2 near the horizon, and therefore that /(r^) = 0. 
Hence we can choose F = y/J in ( 1261) and the solution for the projected Killing potentials 
are _ 

^H^V^^l-/)' (27) 

In the next section we make use of these solutions to express the mass in terms of the 
horizon radius r/j. 

5 Note that if one arbitrarily specifies F at each subsequent time in the evolution, then in general the 
evolution will not be along a Killing field. To keep the evolution along the Killing vector, one imposes 
s a b = K a b = 0, which implies a set of differential equations that F must satisfy. 
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5.2 Dependence of mass on horizon radius 



It is straightforward to evaluate the boundary integral for the ADM mass ( 1131) in terms 
of the coefficient c r that characterizes the far field behavior of the spatial metric. With 
the AdS boundary conditions (12 1|). one finds that the terms in (fT3l) that depend on the 
derivative of the lapse F cancel, and the integral becomes 

M= 16, G (li-l) - (28 > 

where we have taken b\ — 1. We next turn to the task of finding the fall-off coefficient c r 
and hence also the mass in terms of the horizon radius by applying the results ( TT6|) and 
( TT7|) . The integrand at large r in (fT6|) becomes 

r p- 2fc -i) [(1 _ f) u _ (1 _ „ (29) 

while the boundary term on the horizon is easily evaluated using the condition /(r^) = 0. 
Assembling these pieces one arrives at the expression for the thermodynamic potentials 

8 W = -7*^- 2 (rf- 2 *- 1 + (-l) fc ^ry ) (30) 



where is the area of a unit D — 2 sphere. For k = this reproduces the result of [19] 

that 0(°) is given by minus an effective volume of the region behind the horizon. Finally, 
requiring that the sum rule ( fTTl) be satisfied yields the relation 



tr) S -sfeii < 31 > 

/ fc=0 r ft 

This is a key result since it relates the far field behavior of the black hole solution to the 
horizon radius, without making use of an exact analytic expression for f(r) in the region 
between. 



There are two important points to be made about the expression for c r . First, we have 
noted that the solutions for static Lovelock black holes are known to be specified in terms 
of solutions to a polynomial equation of order k max in the function f{r) (9j. This equation 
arises from integrating the Hamiltonian constraint and c r is a constant of integration. 
Equation ( I3TT) may also be obtained by evaluating this polynomial at [52J. What we 
have learned from our more general treatment is that the relation ( 13~T|) expresses the fact 
that the Hamiltonian constraint is a total divergence when defined with respect to evolution 
along a Killing field. Put differently, if a metric does not have a Killing field we have no 
reason to expect that the mass is simply a function of data on the horizon, but in general 
will depend on volume integrals as well. 
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Second, we can ask whether there a computational advantage to working with the Loveloock 
potentials and the sum rule (ITTj) compared to simply writing out the field equations and 
analyzing them? In the spherically symmetric case this is likely a matter of taste. The 
analysis in terms of the Killing-Lovelock potentials is more complicated, but exposes an 
underlying geometrical structure. For more complicated spacetimes, such as rotating black 
holes, it may be that using the Killing potentials allows one to find geometrical relations 
of interest more simply than through a brute force analysis of the field equations. 

We can now substitute in for c r using f[3Tj) and obtain the sought after, result for the mass 
of a static Lovelock black hole in terms of its horizon radius 



M - n D-i r h 3 V- b k ( o 2) 

h h 



This generalizes to all Lovelock theories the result for Gauss-Bonnet black holes given in 
[12J. Such a generalization was assumed to hold in [52] and we have now shown that this 
is indeed the case. 



To get an idea of how this formula behaves, let us examine the dependence of the mass 
on r/j in different limiting regimes. The behavior of the mass as — > turns out to be 
interesting in that it differs between even and odd spacetime dimensions. Since the small 
black holes are dominated by the highest curvature terms, behavior in this limit depends 
on whether the order k max term actually appears in the Lagrangian with non-zero coupling 
constant. We will assume that this is the case. Note that this includes Gauss-Bonnet 
gravity in D = 5 and in D = 6, for which k max = 2, but not in higher dimension. 

By a "small" black hole will be one such that the horizon radius is sufficiently small that 
r 2{k max -k) ^ \b kma Jb k \ for all k < k max . After noting that k max = (D - l)/2, and 2k max = 
(D — 2)/2 respectively in odd and even dimensions, one finds that in the small black hole 
regime the mass depends on as 

M ~ w?G D (D-i) bk ™*x r h D even (33) 

M » (**~ + ^™*-i r l) D odd (34) 

However, for D odd the mass goes to a nonzero value, We see that for even dimensions, the 
mass goes smoothly to zero with the horizon radius as it does for Schwarzschild black holes. 
However, for odd dimensions there is a minimum mass for black holes that is proportional 
to bk max . This minimum mass, or mass gap, has been discussed previously for D = 5 Gauss- 
Bonnet black holes in [12j and for Chern-Simons-Lovelock theories, which have a unique 
constant curvature vacuum, in [15]. 



So long as the coefficient b of the cosmological constant term in the Lovelock Lagrangian is 
non-zero, black holes in the opposite regime of large are dominated by the cosmological 
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constant and look qualitatively the same in all dimensions. A "large" large black hole will 
be one such that the horizon radius satisfies rf^ 3> \bk/bo\ f° r all A; > 0. For large black 
holes one finds that the mass depends on the horizon radius as 



5.3 Surface gravity of Lovelock black holes 

For a metric of the general form ([1]) the surface gravity is given by k = (1/2)0' (r/j). In the 
appendix we argue that near the horizon of a black hole the metric functions must satisfy 
g u ~ —l/g rr , with equality for the functions and their first derivatives at the horizon. 
Hence /(r^) = and /'(r^) = 2k. Evaluating the Hamiltonian function in ( 123]) at 
the horizon and applying the sum rule (IT7j) yields a general relation between the surface 
gravity and the horizon radius without the need for the explicit form of the metric function 
/(r), 

k h v 1 ^ (D - 2k - l)b k . . 

<l^^)) = ^T h l^ — -*i=i) ( 36 ) 

fc=0 'h k=0 'h 

For Gauss-Bonnet gravity this agrees with the expressions for the surface gravity obtained 
from the explicit solutions in [121 S3 E2] using an approach similar to the one here. 

Let us examine at the behavior of the surface gravity in various limits regimes starting with 
small black holes as defined above. In this regime, one finds that 

K ~ ttt - — , D even (37) 

k » - bk ^ rh , D odd (38) 

Again, we see a qualitative difference between even and odd dimensions. In even dimen- 
sions the surface gravity diverges in the limit of vanishing horizon radius, as it does for 
Schwarzschild black holes in D = 4. However for generic Lovelock theories in odd dimen- 
sions, i.e. those in which bk max is nonvanishing, the surface gravity goes smoothly to zero 
with the horizon radius. This change in behavior was noted for Gauss Bonnet gravity in 
D = 5, 6 in [12]. Here we note, in agreement with the observations of [52], that this is 
characteristic of Lovelock black holes in general. 

On the other hand, the surface gravity of large black holes has qualitatively the same 
behavior for all AdS-Lovelock black holes. One finds that 

~ b ° rh fon\ 

« ~ —r^ r 39 

2(D - 2) V ; 

where we have set b\ — 1, so that b\ — (D — 1)(D — 2). In any dimension, therefore, there 
exists a large black hole solution for sufficiently high temperature. In even dimensions also 
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always exists a small black hole at suficiently high temperatures. Hence for D even, and so 
long as the Lovelock couplings are such that the surface gravity stays positive for all horizon 
radii, there will be a minimum temperature at which static black holes exist, while in odd 
dimensions there is no minimum temperature. In addition, there may be local extrema of 
the temperature depending on the choices of the bk- 



5.4 Entropy, thermodynamic potential and Smarr relation 

The expression for the entropy [12] given in section (j5J) is a sum of integrals of Lovelock 
invariants constructed from the induced metric 7 a f, on the horizon cross section. For a 
spherically symmetric black hole 7^ is the metric of a round D — 2 dimensional sphere of 
radius and the Riemann tensor is given simply by R a b cd — (l/ r D^ab- Evaluating the 
various Lovelock terms explicitly gives C^'^^ab) = mE^jk r h 2 ^ k ~^ > - The entropy for a 
static, spherically symmetric Lovelock black hole can then be written as 

swr^- th ( 4o) 



We can also compute the quantity S' defined in (120]) which appears in the Smarr formula. 
It is given in terms of the horizon radius by 

, ^Wf- 2 2fc(fc-l)S fc 



Our goal in this section has been to develop expressions for the thermodynamic properties 
of AdS-Lovelock black holes purely in terms of the horizon radius. The final element of the 
Smarr relation ( fl8|) is the overall thermodynamic potential G, which we will now compute. 
With this in hand, we can check all of our results are consistent with the Smarr relation. 
Substituting our result (l3T|) for c r into the expression (!30~j) for the 9^ and performing the 
sum over k in ( 1201) gives the result 

8^(2, _ 1)2^* 1+ s ( J r f-) (42) 

Combining the results for the mass (132|) . surface gravity (!36|) . entropy (|40l) . S' in (|4T|) and 
thermodynamic potential f)42p . it is now straightforward to check that the validity of the 
Smarr formula (1181) for the AdS-Lovelock black holes. 
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6 Free energy and phase transitions 



In this section, we will make use of our results to give a general expression for the free 
energy of AdS-Lovelock black holes. We will discuss the behavior of this expression in the 
limits of large and small horizon radius, in which it is independent of the detailed choice 
of Lovelock couplings. 

The free energy of asymptotically AdS black holes has been a topic of interest since the 
work of Hawking and Page [56] in Einstein gravity. They found that there exists a minimum 
temperature To for such black holes which occurs for a horizon radius tq. The temperature 
diverges both in the limit of large black holes and in the limit of zero horizon radius. Black 
holes with horizon < r , like asymptotically flat black holes, have negative specific heat 
and cannot be in stable equilibrium with a thermal bath of radiation. However, solutions 
with rh > r have positive specific heat and can be in stable equilibrium. For large black 
holes with temperatures just above T , the free energy is positive and thermal AdS space, 
with zero free energy, represents that globally preferred thermodynamic state. However, 
the free energy of large black holes becomes negative above a critical temperature T\ > T 
(or correspondingly for black holes with radii exceeding a certain threshold r\) and the 
black hole is then the globally preferred state. 

Hawking and Page found the free energy by computing the Euclidean action, which re- 
quires the analytic form of the AdS-Schwarzschild spacetimes. Similar computations of 
the free energy have been carried out in Gauss-Bonnet gravity [121 SSI EZl EE] , where the 
explicit solutions are also known [TUl [EE] . An expresssion for free energy in general Lovelock 
gravity is given in [59J which is calculated using a generalized quasi-local mass defined for 
spherically (plane or hyperbolic) symmetric spacetimes. However, theses computations are 
considerably more complicated. Our method offer a much a simpler route to the result 
and offers some degree of physical interpretation to the different terms in the result. It 
also yields the answer in the general Lovelock case, where do to the absence of explicit 
analytic solutions, computation of the Euclidean action may not be practical. Moreover as 
stated above, we envision further applications to even less well understood solutions, such 
as rotating Lovelock black holes. 

It is worth noting explicitly that the equality of the free energy with the Euclidean action 
times the temperature continues to hold generally in Lovelock gravity theories. Hawking 
and Horowitz [HD] used the Hamiltonian framework for Einstein gravity to demonstrate that 
Ie = f3M — S where /3 is the Euclidean period. A similar construction works in Lovelock 
gravity, the basic steps being as follows. 

Write the volume term of the Euclidean action for a static black hole in Hamiltonian vari- 
ables, and then directly derive the boundary term for the action by varying this expression. 
The variation of the volume term is the same calculation that is done to find the mass (see 
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[T5j[20]) with an additional integration over Euclidean time. Evaluated on static solutions, 
the volume term of the action reduces to a sum of the constraints and hence vanishes. The 
value of the action is then given by the boundary terms. At infinity the boundary term is 
simply M/T. The Euclidean metric does not have a horizon. However, Hamiltonian evolu- 
tion with respect to Euclidean time fails to be well defined at r^. In order to compensate 
for this, one introduces an inner boundary at r = + e. In the limit e -)• the boundary 
term is equal to the entropy S. One then has the result that the free energy defined by 
F = TIe coincides with the thermodynamic free energy F = M — TS. Note also that the 
free energy has been defined such that it vanishes for 

Let us see how the Hawking-Page transition arises in the present framework. In Lovelock 
gravity the cosmological constant is given by A = —bo/2, so that 60 is positive for A 
negative. Starting from F = M — TS, let us use the Smarr formula fll8p to eliminate the 
mass. For Einstein gravity the quantities sp) and S' = in the Smarr relation vanish, 
while as shown in (19] the overall Lovelock thermodynamic potential is given by = —Vbh 
where Vbh is an effective volume for the black hole given by the flat (or AdS) spacetime 
volume of a sphere of radius r^. One then arrives at an expression for the free energy 



The important observation is that the apparent contribution of kA to the free energy has 
changed sign in (|4*3|) relative to the original expression F = M — TS. This is because the 
mass receives a positive contribution from kA in the Smarr formula and the overall net 
coefficient is always positive. For an asymptotically flat black hole, i.e. with vanishing bo, 
this is the only term and the free energy is always positive. With a negative cosmological 
constant there is a negative definite contribution that takes the form of a cosmological 
pressure bo times an effective volume of the black hole Vbh- Hence the Hawking-Page phase 
transition at which the free energy changes sign may be thought of as arising from the O 
term in the Smarr relation, which is itself analogous to a P\^-type contribution in classical 
thermodynamics. 

To determine whether the positive or negative term dominates F one needs to use know how 
the surface gravity behaves. Hawking and Page used the analytic solutions to compute k 
and found that the free energy is positive for small black holes as in Schwarzchild. However, 
they found that large black holes have negative free energy. Indeed, substituting k from 
equation (I36p one recovers these results. However, rather than recall this case in more 
detail, we turn to an analysis of the free energy for general Lovelock black holes, which 
includes the case of Einstein gravity. We will see that the behavior of the free energy for 
large black holes is qualitatively the same as the Hawking-Page case. On the other hand, 

6 However, because of the presence of a mass gap for odd dimensional AdS-Lovelock black holes, the 
limit Th — > will not vanish in our subsequent expressions for F given below. 



1 



(kA - boVbh) 



(43) 



F 



8nG(D - 3) 




(44) 
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the free energy of small black holes differs between even and odd dimensions due to the 
alternating behavior of the surface gravity. 



6.1 Thermodynamic stability and phase transitions for Lovelock 
black holes 



There are a number of contributions to the general Lovelock free energy and the simplifi- 
cation made above in the Einstein case using the Smarr formula no longer yields an easily 
interpretable expression. Instead we simply substitute into the free energy the formulas for 
the mass (|32|) and the entropy (HUT) to obtain 



^D- 2 rf- 3 ^ b k / 2k 



- iQnG( D - 1) £ ^ y 1 Krh w^k)) (45) 

where k can also be regarded as in ( )36|) as a function of the horizon radius and the Lovelock 
couplings. This expression for F agrees with that of reference [12j in the case of Gauss- 
Bonnet gravity and [5H] in general Lovelock gravity, but disagrees^ with that given in 
reference 



Let us start by examining the behavior of the free energy in the large black hole limit. 
Form fl39|) . we see that the surface gravity grows like &o r /i for large in all dimensions. 
One then finds that there are both positive and negative contributions at leading order 
growing like b^r®' 1 . The net result for the free energy in the large black hole limit turns 



out to be negative, 



° 2 ° h (46) 



!6irG(D- l)(D-2) 

This result was also pointed out in [61] (see also [621 163] ) and is not surprising since the 
behavior in this regime is dominated by the cosmological constant and Einstein terms. 

The behavior of the free energy ( 1451) in the small black hole regime, on the other hand, is 
dominated by the highest curvature terms and differs between even and odd dimensions. 
In even dimensions both the mass M and the product kS scale like in the limit of small 
horizon radius. One finds in this case that the positive contribution to the free energy 
coming from the mass dominates, giving 

„ fin_o bu Th „ . . 

W 32irG(D-l) ' D 6Ven (4?) 



7 The expression for the free energy in [13] is obtained by computing the volume term in the Euclidean 
action, apparently without the inclusion of a boundary term. A subtraction of the action for pure AdS 
at large radius is used to regularize the result. There are then several differences with our calculation. 
First, we implicitly use the Euclidean action with the boundary term that gives a well defined variational 
principal in the Hamiltonian variables. As described above the action on solutions is then given entirely 
by the boundary term and is equal to the finite quanity M — TS, without a need for regularization. 
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In odd dimensions both the surface gravity and entropy vanish as — > (see equation 
(155]) ). while the mass has a finite positive limiting value, giving 



F « n °- 2 ^ r/ V D odd (48) 

Hence F is positive for small black holes in all dimensions. For D even F goes to zero, 
while for D odd F goes to a nonzero positive value. 

It is also straightforward to compute the specific heat in the large and small black hole 
limits. Using the expressions for the mass and the surface gravity in section (JSJ) we find for 
large black holes 

8M (D-2)n D ^ 2 r°~ 2 



dT AG 



(49) 



and for small black holes 



9M kmax n D - 2 b kmax r 2 h 



_ Deven (50) 

w " Dodd (51) 

where we have assumed that the couplings bn, hh and hh _i are all nonzero. 

sr o yJi Fornax ^max - 1 - 

Let us now summarize these results. While a detailed understanding of the global and 
local thermodynamic stability of AdS-Lovelock black holes throughout the entire range of 
horizon radii would require specifying the entire set of Lovelock couplings, the behavior in 
the large and small black hole regimes is generic. In dimensions, we find that the behavior 
of AdS-Lovelock black holes in these regimes is similar to that found by Hawking and Page 
in Einstein gravity [56J . Black holes become arbitrarily hot in both limits. Small black holes 
in even dimensions exhibit negative specific heat and have positive free energy, indicating 
instability to both perturbative and non-perturbative fluctuations. Large black holes, on 
the other hand, have positive specific heat and negative free energy, both indicating their 
stability to small thermal fluctuations and that they are the thermodynamically globally 
preferred state. If one assumes that the Lovelock couplings are such that the temperature 
stays positive for all r^, then there must exist a minimum temperature T below which no 
black hole solutions exist. As a consequence of these similarities, we can expect that even 
dimensional Lovelock black holes will at least have a simple Hawking-Page phase transition 
and possibly a more complicated structure of phase transitions, depending on the detailed 
behavior of the temperature and free energy over the whole range of horizon radii. 

In odd dimensions only large black holes exist at very high temperatures. They have 
positive specific heat and negative free energy as in even dimensions and are therefore 
thermodynamically stable. The odd dimensional small black holes, on the other hand, 
have positive positive free energy and also have positive specific heat. So they are stable to 
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small, but not large, thermal fluctuations, and unlike small black holes is even dimensions 
can be in stable equilibrium with a thermal bath. 

In odd dimensions the low temperature picture is different. This was studied in D = 5 
Gauss-Bonnet gravity in references [T2], [58] . Black holes exist for [12] temperatures down 
to T = 0, and the low temperature black holes have positive specific heat, so there is a 
locally stable alternative to the gas state. At these low temperatures the pure gas state 
has lower free energy since F gas — > as T — > 0, while the black hole starts with F of order 
bk max - Still, a small black hole can exist in equilibrium with a low temperature gas, unlike 
the situation in even dimensions. So the behavior of small black holes in D = 5 Gauss- 
Bonnet gravity continues in odd dimensions as long as the highest curvature Lovelock term 
is included. For very high temperatures there is only the one large black hole state. As in 
even dimensions this is both globally and locally preferred. 



7 Discussion 

In this paper we started by deriving a fully explicit formula for the ADM mass of an 
asymptotically AdS spacetime in a generic Lovelock gravity theory, via the Hamiltonian 
methods of Regge & Teitelboim [2]. We then proceeded to study various thermodynamic 
properties of AdS-Lovelock black holes. In particular, we made use of the Killing-Lovelock 
potentials that exist in these spacetimes in order to evaluate the mass in terms of the horizon 
radius and the Lovelock couplings. After finding expressions for the surface gravity and 
entropy, these ingredients were assembled to give the free energy. All of these expressions 
are quite general, assuming only that solutions exist with the prescribed asymptotic forms. 

As mentioned in the introduction, we envision further applications of these techniques to 
stationary solutions that are even "more unknown", such as higher dimensional rotating 
charged black holes, or rotating Lovelock black holes. Another possible direction for future 
work would be to look at black holes/branes in AdS with planar horizons. In this case 
the asymptotic boundary of a spatial slice is a plane rather than a sphere. If one of these 
directions is compact with length L, then it would be necessary to further extend the first 
law to include an appropriate 5L term as in the asymptotically flat Kaluza-Klein case [61] . 
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Appendix - Near horizon behavior of g rr 

The horizon of a static black hole occurs where the norm of the static Killing field is zero. 
In the coordinates of the metric ^ this is at such that 4>(rh) = 0. In this appendix 
we show that if one assumes that the surface gravity is finite and that the Hamiltonian 
constraint equation is satisfied at r^, then near the horizon it follows that f(r) ~ <j){r), and 
in particular that /(r^) = , f'frh) = 2/t. 

We point out that for the spherically symmetric black hole metric ([1]) it has been shown 
[IT] [9] that the vacuum field equations imply that —g tt = ^/g rr , or 0(r) = f(r) every- 
where. Clearly this is a stronger result than the result shown here. However, to derive the 
expressions for M, k, and the 0^ we only need to use the values of / and /' at r^, which 
can be found with correspondingly less work, as follows. 

Assume that the metric function in (pQ) goes to zero like a power law as r approaches r^, 
i.e. that near the horizon 

<j>^<h.{r-r h Y (Al) 
where (j>x is a constant. The surface gravity, given by k 2 = — |(V a ^)V a ^ b , then becomes 

k 2 = -p 2 lim (r - r H ) p ~ 2 f{r) (A2) 

2 r->r h 

In order for k 2 to be finite, it must be that the metric function f(r) behaves like 

/~A(r-r h ) 2 - p (A3) 

with fx another constant. Let us rewrite the expression PU for H ( l ] as 

H (k) = r i-2k [ (jD _ 2k _ _ /)/r _ fc(1 _ (A4) 

and consider the constraint equation f !24p . If the power law index p > 1, we see from line 
( ]A3j) that /' diverges at the horizon, and inspection shows that H = cannot be satisfied 
at rh- On the other hand, if the power law index p < 1, then /(rv,) = f'(rh) = and again 
the constraint cannot be satisfied^] at r^. Hence, in order to have finite, non-zero surface 



8 If the spacetime is not vaccuum, the the gravitational constraint is H = — Wirp, and so powers p < 1 
are not ruled out. 
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gravity, the power law index in (1A2|) and (1A3j) must be p = 1, so that ~ 0i(r — 77J and 
/ ~ /i( r — r /i); which gives k 2 = We see that the overall numerical scale of k is not 

fixed by this argument, which makes sense because the normalization of k is fixed by the 
norm of the time-translation Killing vector at infinity, which requires knowing the function 
0(r) throughout the spacetime. However, one can fix the scale by requiring that it gives 
the right answer for Schwarzchild. 
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